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Abstract 

We study the universality between a discrete spin model with icosa- 
hedral symmetry and the 0(3) model in two dimensions. For this 
purpose we study numerically the renormalized two-point functions of 
the spin field and the four point coupling constant. We find that those 
quantities seem to have the same continuum limits in the two models. 
This has far reaching consequences, because the icosahedron model is 
not asymptotically free in the sense that the coupling constant pro- 
posed by Liischer, Weisz and Wolff does not approach zero in the 
short distance limit. By universality this then also applies to the 0(3) 
model, contrary to the predictions of perturbation theory. 
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1 Introduction 



The subject of the enhancement of a discrete symmetry to a continuous one 
at large distances is important both theoretically and phenomenologically. 
Indeed, even if there are good grounds to expect that a certain material 
is well described by some model enjoying O(N) symmetry, one may wonder 
what might be the effect of anisotropies. This question was addressed in 1977 
by Jose et al p[ for the 0(2) nonlinear a model. Nonrigorous renormalization 
group arguments led them to the conclusion that in two dimensions (2D) the 
discrete symmetry should be enhanced to full 0(2) invariance if N > 5 for 
f3 (inverse temperature) not too large. For N sufficiently large, the occurrence 
of this phenomenon was proven rigorously by Frohlich and Spencer in 1981 
0, who showed that there exists a range of temperatures in which spin 
correlation functions decay algebraically and are 0(2) invariant. Frohlich and 
Spencer proved also that a similar phenomenon of discrete abelian symetry 
enhancement occurs in 4D gauge theories. 

For a long time, the consensus was that no symmetry enhancement should 
occur in nonabelian models. The main reason appears to have been the be- 
lief that for continous symmetries, these models exhibit asymptotic freedom 
(AF). The discrete models, known rigorously to undergo phase transitions at 
nonzero temperature, did not seem likely to be AF, hence had to be differ- 
ent. A proposal for nonabelian symmetry enhancement came however from 
Newman and Schulman [|J . Their argument was based on the fact that if the 
discrete symmetry group is sufficiently large, any 4-th order polynomial in- 
variant under the discrete group is also invariant under the continuous group 
in which it is contained. While this is an undisputable mathematical fact, the 
question was why 4-th order? Their heuristic answer was that the renormal- 
ization group flow was expected to be free of bifurcations as the dimension 
D was varied between 2 and 4. So if one started just below D — 4, where 
the most one could have is a 4 interaction (higher powers being irrelevant), 
this symmetry enhancement should persist down to D — 2. Since however 
in D = 2 all polynomials in <ft are relevant and it is easy to write down 
such polynomials possessing only a discrete symmetry, and to construct the 
corresponding P(<p)2 models, the validity of the argument of Newman and 
Schulman remains unclear. 

A different heuristic argument in favor of symmetry enhancement for 
abelian as well as nonabelian groups was put forward by Patrascioiu in 1985 
[J5J. For spin models, his argument went as follows: at sufficiently low temper- 
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atures, there exists a phase with long range order (l.r.o.) because, as Peierls 
showed long ago, given an ordered state, there is not enough free energy to 
create a domain in which the spin points elsewhere. Now consider a model 
like Z5. As one increases the temperature, clearly the first abundant domains 
to form would be those in which the spin pointed in a direction immediately 
neighboring the one chosen by the boundary conditions (b.c.) for the ordered 
state. For temperatures not too high, the system could form domains inside 
domains of neighboring spin values. This would be different from the phase 
at high temperature, where no such restriction between adjacent domains 
would be required. This scenario does not seem to have anything it do with 
the model being abelian or not and Patrascioiu suggested that, since it was 
known to happen in abelian models, it must happen also in nonabelian cases. 

Except for the papers quoted above, in the 80s, while it was quite fash- 
ionable to replace continous groups with discrete ones in Monte Carlo simu- 
lations, everybody seemed to be convinced that the discrete and continuous 
models belonged to different universality classes. Our interest in the subject 
was rekindled in 1990 when, together with Richard, we derived a rigorous 
inequality relating correlation functions in the dodecahedron model to those 
of Zio P|. Our result was that for any (3 the dodecahedron model is more 
ordered than Z\q at 0.60 7 2 /?. Since it is was pretty well established that 
Z w possess an extended intermediate phase which is 0(2) invariant, our 
inequality implied that provided /3 m (D) > (3 C (Z W ) the dodecahedron must 
also possess an intermediate massless phase. Here (5 m (D) denotes the onset 
of the l.r.o. phase in the dodecahedron and j3 c (Z w ) the onset of algebraic 
decay in Z w . We determined numerically these values @ and concluded 
that the dodecahedron seemed to possess an intermediate massless phase 
for 2.15 < (3 < 2.8. We conjectured that this phase must enjoy full 0(3) 
invariance. 

Intrigued by our findings for the dodecahedron, in the early 90s we looked 
numerically at the other regular polyhedra. While the cube is obviously 
equivalent to 3 uncoupled Ising models, hence not a good candidate for ex- 
hibiting 0(3) invariance, the other 3 regular polyhedra (platonic solids) a 
priori are. Actually in the scenario advocated by Patrascioiu in 1985 ||, the 
tetrahedron should not be able to simulate spin waves since its spin gradient 
can take only one nontrivial value (in fact it is nothing but the 4 state Potts 
model). The octahedron and the icosahedron could. Our numerics suggested 
that the octahedron had a first order transition. For the icosahedron the MC 
data suggested a second order transition from the high temperature phase to 
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the low temperature phase exihibiting l.r.o.; in particular there did not seem 
to be an extended massless phase, as for the dodecahedron. 

That the discrete icosahedral symmetry may be enhanced to 0(3) began 
to become manifest in 1998 when we began extensive numerical investigations 
of the contiuum limit of the spin 2-point function versus p/m || in the 
dodecahedron model. The original motivation of that study was to compare 
the lattice continuum limit with the form factor (FF) prediction of Balog and 
Niedermaier . Since it was not to be expected that the latter could possibly 
describe the continuum limit of the dodecahedron model, we decided to take 
data on this model too. To our surprise, the 0(3) data seemed to agree 
with both the FF prediction and with the dodecahedron. Since Balog and 
Niedermaier had produced a convincing argument |TD[ that the FF approach 
incorporates AF and we could not see how a discrete model, freezing at 
nonzero temperature, could possibly be AF, we decided to refine our data 



by concentrating on the region p/m < 13. Our results ]TT[ showed small 
but statistically significant deviations between 0(3) and FF, but excellent 
agreement between 0(3) and the dodecahedron. 

The comparison of the FF and 0(3) could have been marred by lattice 
artefacts, a problem to which we will return below. This is why a different 
comparison was performed by Patrascioiu |l2j , who computed the renormal- 
ized spin 2-point function versus the physical distance x/£. The results were 
very similar with the ones we will report here about the icosahedron. They 
suggested that the dodecahedron and 0(3) models share the same continuum 
limit. 

Another reason to investigate the discrete spin models arose in 1999, 
while in collaboration with Balog, Niedermaier, Niedermayer and Weisz we 
decided to compare the FF prediction for the renormalized coupling (to 
be defined below) with its lattice continuum limit value. Although our first 
results |13[ suggested excellent agreement, as we continued to reduce the 
error bars and especially to take data at larger correlation length £ rs 167, 
our original extrapolation to the continuum limit became dubious and in our 
second paper |TJ| we stated that we could no longer give a reliable number 
for the lattice continuum limit value of g^. In the hope of gaining insight 
into this issue, we turned again to the discrete models. Corroborating our 
previous findings regarding G r (p/m) and G r (x/£), the data, to be shown 
later, suggested that qr also agreed in the dodecahedron, icosahedron and 
0(3) models. 

In the course of the debates of our collaboration F. Niedermayer raised 
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the intriguing possibility that maybe indeed all these models do have the 
same continuum limit, which however is AF. The present paper is our reply 
to his suggestion. The results were communicated to him already in 1999, 
hence we were surprised by his recent paper with P. Hasenfratz ||15|| . In 



their paper, Hasenfratz and Niedermayer claim that indeed they find strong 
numerical evidence that the dodecahedron and icosahedron models have the 
same continuum limit as 0(3), but state that, since "overwhelming evidence 
exists that the 0(3) model is AF", the dodecahedron and icosahedron models 
must be AF too. The authors do not mention which 'overwhelming evidence' 
they have in mind, but omit the following facts, which in our opinion are 
relevant to this issue: 



The data for gn in 0(3) produced by Balog et al at larger £ values 
do not support the l/£ fit for the lattice artefacts. 



The last report of Balog et al [|14| retracted the original prediction 



gn = 6.77(2) and stated that due to uncontrollable lattice artefacts, no 
continuum value could be given. 

• In our recent paper |16[], we combined mathematically rigorous argu- 
ments with some numerics to conclude that the 0(3) model must un- 
dergo a transition to a massless phase at finite (3. We then proved 
rigorously that such a phase transition rules out the existence of AF in 
the massive continuum limit. 

Being based on numerics, any of the above stated facts could be false, but if 
Hasenfratz and Niedermayer have evidence to that effect, they should put it 
forth so that it can be scrutinized. 

In this paper we will compare the continuum limit of the 0(3) model 
to that of the icosahedron model by comparing the values of g^ and of the 
renormalized spin 2-point function G r (x/^). The advantage of the icosahe- 
dron model is the existence of a rather well localisable critical point, whereas, 
as stated above, the dodecahedral model appears to have a soft intermediate 
phase. To address the issue of AF, we study the value of the Liischer-Weisz- 
Wolff (LWW) coupling constant at the critical point. According to LWW, 
AF requires that the continuum limit of this observable vanishes at short 
distances. If in fact the icosahedron model has the same continuum limit as 
0(3), then the LWW coupling constant should vanish at the critical point 
of the former model. We find excellent evidence that it does not. Thus our 
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conclusion is that either, in spite of the excellent agreement observed for 
£ ~ 121, 0(3) and the icosahedron model have different continuum limits, or 
neither is AF. 

The paper is organized as follows: we first describe the critical properties 
of the icosahedon model and in particular locate its critical point. This allows 
us to determine the value of the LWW running coupling at the critical point. 
Then we move on to the comparison of the icosahedron and the 0(3) models. 
We compare the renormalized coupling constant and the renormalized spin- 
spin correlation function of the two models and present convincing evidence 
that they converge to the same continuum limit. 



2 Critical Behavior of Icosahedron Model 

We first describe the critical properties of the icosahedron model: we locate 
the critical point and give some estimates of critical exponents. This allows us 
to determine the value of the Luscher-Weisz- Wolff (LWW) coupling constant 
at the critical point, which is independent ot the size L of the lattice, as 
required by scaling. This (nonzero) value is also the short distance limit of 
the continuum value of the LWW coupling constant. 

The icosahedron model is defined by the standard nearest neighbor cou- 
pling between the spins 

H=-Y,8i-8j (1) 

(ij) 

where the spins s$ are unit vectors forming the vertices of a regular icosahe- 
dron. In suitable coordinates, those 12 vertices are given by 

Ijtk 2irk 

e k = (scos^-,ssin^-,ccos(7rA;)) (k — 1, 2, ...10); (2) 

where 
and 

en = (0,0,1) and e 12 = (0,0,-1) (4) 

The model has at least two phases, a high temperature phase with ex- 
ponential clustering and full symmetry under the icosahedral group Y, and 
a low temperature phase with spontaneous magnetization and 12 coexist- 
ing phases, with the magnetization pointing into one of the 12 directions of 
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the icosahedron. At intermediate temperatures there could be in principle a 
phase with partial breaking of Y, but there is no reasonable candidate for 
a possible unbroken subgroup. There is also the possibility of an extended 
intermediate phase with no symmetry breaking, but actual enhancement of 
the symmetry from Y to 0(3) as well as only algebraic decay of correlations 
analogous to the Z N models and the dodecahedron model (see above). 

In the icosahedron model the situation seems to be simpler, however: it 
seems to have a single critical point separating the high and low temperature 
phases, as we will demonstrate. 

To determine the critical point we proceed as follows: we determine the 
mass gap m(L) = 1/£(L) in an (ideally infinite) strip of width L; in practice 
we use a finite strip of size L x L t with L t » L and measure the effective 
correlation length defined by 



where 




We then study the behavior of the quantity 

- - 2L 

9 = (N-mLy (8) 

which we consider as a function g(z,£) of z = L/£(oo) and the infinite volume 
correlation length £(oo). The continuum limit of this quantity for fixed z is 
the 'running coupling constant' introduced by Luscher, Weisz and Wolff M 
for the 0(N) models. In the high temperature phase, where the model has a 
mass gap m(oo) = l/£(oo) in the infinite volume limit, g will grow linearly 
with L. On the other hand in the low temperature magnetized phase, the 
mass gap in a finite volume goes to faster than 1/L, so g will decrease to 
as L —>■ oo. At a critical point the model will be scale invariant for large 
distances and g should converge to a finite nonzero limit. 

To determine g we took data on lattices of size L x L t with L t = 10L; L 
was varied from 10 to 320. Fig.|l] (Table 1) shows clearly the dramatic change 
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Figure 1: The LWW coupling constant g as a function of L for various values 
ofP 



in behavior between (3 = 1.800 and (3 = 1.805. For (3 = 1.802 g shows only 
some small variation for small L (L < 40) and stabilizes for larger L. So we 
estimate 

P crt = 1.802(1) (9) 

where the error is of course somewhat subjective. 

To corroborate this determination of the critical point, we also measured 
on the same lattices the renormalized coupling defined as 

/5 g 4 \ LL t 

Here g^ is the magnetic susceptibility multiplied by the volume of the lattice 
L x L t , i.e. 

92 = Y,( s i ■ s i) ( n ) 

and 

#4= (( s h ■ s fa)( s i 3 • S *J> ( 12 ) 
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Figure 2: The renormalized coupling for asymmetric L x L t lattices as a 
function of L for various values of (3 

This quantity is also a Renormalization Group (RG) invariant and therefore 
should also go to a constant for (3 = (3 cr t. It is well known that in the high 
temperature phase Qr(L) goes to a nontrivial constant for L —>■ oo, (which 
as about 6.7, see below) whereas in the magnetized phase it goes to 0, so 
one expects a qualitatively similar picture as for g. The data presented in 
Table 1 and displayed in Fig.|2] confirm this nicely and are consistent with 
the estimate of (3 crt given above. 

Finally we want to see if our determination of (3 cr t is consistent with a 
singularity in the thermodynamic values of the correlation length £ and the 
susceptibility x- We therefore measured £ and x on lattices with L/£ ~ 7 
for various values of /3 < 1.802; our data are given in Tab. 2 . There is a 
row listing the number of runs; a run consists of 100,000 cluster updates 
for thermalization followed by 20,000 sweeps of the lattice for measurements. 
Each run is started independently with a randomly chosen new configuration. 

To describe the critical behavior of the data for £ and x two types of fits 
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were tried: firstly a Kosterlitz type fit with an exponential singularity 

^MteM- *^ exp hcN' (13) 

and secondly a power law fit of the type 

£ = Q(&rt -/*)-•', X = C x ([3 C rt-{3)-" (14) 

or similar ones in 1/(3. Both types of fit are not very good, with similar 
quite large values of x 2 an d they do not allow a very precise determination 
of the fit parameters. The reason seems to be the following: the asymptotic 
singular behavior seems to have significant subleading contributions, which 
however cannot be well determined with only 5 values of (3. Trying to fit 
our very precise data with functions that do not describe the behavior with 
similar accuracy necessarily leads to a poor fit quality, even though visually 
the data may be very well described (see Fig.^j) . 

The Kosterlitz type fit leads to predictions of (3 crt ~ 1.93, which is unac- 
ceptably large - this value is deeply in the magnetized phase. The power law 
fits, on the other hand, give values of (3 crt quite close to our preferred value 
1.802 . By playing with the number of parameters and fitting in 1//3 as well 
as (3 we obtain quite spread in the values of the exponents; they fall into the 
intervals 

^ = 1.6 - 2.0 and 7 = 2.9-3.5 (15) 
This is consistent with a value of 

7] = 2 - - « .25 
v 

which is also favored by the data for spin-spin correlation function (see be- 
low). 

In Fig.^ we use the best values produced by the power law fit in 1/(3 with 
no subleading corrections, namely 

^ = 1.717, 7 = 3.002. (17) 

We plot ^~ 1 / v and x -1 ^ 7 vs P together with the fits, which are straight lines 
intersecting the abscissa at (3 — 1.802. So even though the thermodynamic 
data for £ and x do not lead to a precise prediction of the critical point and 
the critical exponents, they are certainly consistent with our determination 
based on the LWW coupling constant. 

10 



(16) 



We also investigated the possibility that the transition from the high 
temperature phase to the one with long range order is first order, but we did 
not find any signal for phase coexistence. 




0.55 0.6 0.65 0.7 



Figure 3: Thermodynamic values of £ x l v and \ l ^ with v = 1.717 and 
7 = 3.002 as functions of P; the lines are fits. 

3 The renormalized coupling in the icosahe- 
dron and 0(3) models 

To check whether the icosahedron and the 0(3) model define the same con- 
tinuum limit, we also determined the renormalized coupling constant on 
thermodynamic lattices. For this purpose we took data on square lattices 
with z = L/£ ~ 7. To correct for the small finite size effects still present, 
we determined a finite size scaling curve by taking data at various values of 
z = L/£ at (5 = 1.550 corresponding to £ ~ 19.7. The finite size effects are 
fitted with the function 

9r(z) = g R (oo)(l - dy/zexp(-z). (18) 
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Figure 4: Finite size scaling of qr 



This behavior is inspired by the spherical model and was successfully used 
, rj|. In Fig.[|we show the data at (3 = 1.55 together with the 
fit, which has a \ 2 P er d.o.f. of 1.1 ; the constant d is determined to be 



already in 



d = 4.63 ±.15 



(19) 



The results of this fit are then used to extrapolate our data taken with z « 7 
to the thermodynamic values. In Tab. 2 we present our data together with 
the extrapolation. 

Finally we compare our data for g^ in the icosahedron model with those 
we had obtained earlier for the 0(3) model [1^] and again presented in Tab. 2. 
Fig. [5] shows the thermodynamically extrapolated data for the two models 
for various values of £. Even though the lattice artefacts are quite different 
for the two models, and in spite of the fact that we are not quite sure how 
one should extrapolate to the continuum limit, the data show that the two 
models approach each other with increasing £ (decreasing lattice spacing) 
and suggest that they will have the same continuum limit. 
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Figure 5: Comparison of g R in the icosahedron and 0(3) models Full sym- 
bols:icosahedron; open symbols: 0(3) 

4 Spin correlation function in the icosahe- 
dron and 0(3) models 

In this section we compare the renormalized spin-spin correlation functions 
of the icosahedron and the 0(3) models. They are defined as 

G(i/0 = ^M()) •*(<)> (20) 

The physical distance is 

x= l - (21) 

We measured the two-point functions in the icosoahedron model at (3 — 1.47, 
1.55, 1.61, 1.665 and 1.707 corresponding to £ « 11, 20, 34, 64, 122 (see 
Tab. 2). For the 0(3) model we used the data from ref.|]12[ at (3 = 1.8, 1.9, 



1.95 corresponding to £ ~ 65, 122 and 168. In Fig.|6| we show G(x)x 1 ^ 4 for 
the 0(3) model together with the 2-loop PT prediction M 
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Figure 6: The renormalized spin-spin correlation for the 0(3) model for var- 
ious lattice spacings 



_ 1.000034657 
^ ~ 3vr 3 1.001687 2 



1, .1159 
t + lnt + 1.1159 + -\nt + 



t t j < 22 > 

where t = — ln(x/8) — 1. It can be seen that the data approach their con- 
tinuum limit from above and deviate considerably from the PT prediction. 
For x > 0.4, barring some very slow convergence to the continuum limit, our 
data suggest that the lattice artefacts are quite small for the large correlation 
lengths we are using. 

In Fig.^ we present G{x)x 1 ^ for the icosahedron model together with the 
same expression for 0(3) at approximately the largest value of £ ~ 122. The 
lattice artefacts have the opposite sign at least for the lattices with £ < 122, 
i.e. the data are increasing with decreasing lattice spacing. At £ ~ 122 they 
are already quite close to the corresponding correlation function in the 0(3) 
model. It is therefore reasonable to expect that if we could further refine the 
lattice we would see that the data follow the behavior of the 0(3) model, 
i.e. the lattice artefacts would change sign and the final approach to the 
continuum limit would be from above. 

It should be stressed that already at £ m 122 the icosahedron and 0(3) 
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Figure 7: The renormalized spin-spin correlation for the icosahedron model 
for various lattice spacings: solid line 0(3) model for £ « 122 

data are much closer to each other than the 0(3) data are to the PT predic- 
tion. 



5 Conclusion: Universality between the icosa- 
hedron and 0(3) models and asymptotic free- 
dom 

We have accumulated strong evidence that the continuum limits of the dis- 
crete icosahedron model and the continuous classical Heisenberg (0(3)) model 
describe the same quantum field theory. 

As discussed in the introduction, this is one more fact which puts the 
asymptotic freedom of the 0(3) model severely into doubt. The point of 
view advocated by Hasenfratz and Niedermayer [|15|], namely that the con- 



tinuum limit of the discrete icosahedron model should be asymptotically free 
is untenable in view of our results about the LWW running coupling g: our 
data (see Fig.0) indicate that g(L) runs to a fixed point value g* ~ .59 at 
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small distances. Actually to determine the true running coupling, one should 
take first the continuum limit at fixed z = L/£ and then the limit z — > 0. 
Since this is not feasible, we instead studied the finite size scaling at and 
around the critical point and took as our estimate of g* the apparent limit 

lim m(L)L (23) 

L^oo 

at the critical point. Somebody might wonder if it is not possible that the 
continuum limit shows asymptotic freedom after all; in other words, if it is 
possible that 

lire L lim g(z,£) = ? (24) 

2— »0 £^oc 

Our data make this extremely unlikely: for fixed (3 < (3 crt ~ 1.802 g is 
increasing with L (except for very small lattices) and it is always larger 
than .59, even at (3 = 1.802, our estimated critical point. To claim that in 
the continuum limit g(z, oo) would go to one would have to assume some 
truly bizarre L dependence at fixed (3. This is illustrated in our Fig|| which 
shows some of our data for g as a function of z = L/£(oo). We used data 
at (3 = 1.665 and 1.707, where we know the correlation length £(oo) quite 
well, together with the data taken slightly below the estimated critical point, 
wehre we used the fit appearing in Fig.|| to estimate £(oo). The solid curve 
is a fit of the form 

g( z ) = g * + az 1/2 + bz (25) 

Since we did not make any effort to control the lattice artefacts and 
estimate the precise continuum values, this figure should be taken with some 
caution. It does, however, illustrate nicely the qualitative behavior of the 
LWW running coupling near the critical point. 

To sum up: the universality observed between the icosahedron and the 
0(3) model gives strong evidence against asymptotic freedom of the latter. 
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Figure 8: The LWW running coupling for the iscosahedron model as a func- 
tion of L/£(oo) 

Tab.l: 

The LWW running coupling g and the renormalized coupling g^ as a function 
of L for various values of (3 



f3 


= 1.665: 






L 


10 


20 


40 


9 


.9121(14) 


1.0651(11) 


1.3033(21) 


9r 


2.216(13) 


2.501(12) 


3.080(29) 



(3 = 1.707: 



L 


10 


20 


40 


80 


9 


.8156(9) 


.9312(12) 


1.0882(23) 


1.3276(35) 


9r 


1.898(6) 


2.213(11) 


2.593(23) 


3.171(36) 


P 


= 1.75: 








L 


10 


20 






9 


.7221(24) 


.7902(21) 






9r 


1.695(14) 


1.841(24) 
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p 

L 


= 1.80: 

10 


20 


40 


80 


160 




9 
9r 


.6083(16) 
1.363(9) 


.5989(17) 
1.329(9) 


.6018(21) 
1.357(14) 


.6158(17) 
1.395(9) 


.6337(22) 
1.459(13) 




P 

L 


= 1.802: 

10 


20 


40 


80 


160 


320 


9 
9r 


.6015(10) 
1.340(4) 


.5951(14) 
1.3330(6) 


.5919(26) 
1.328(5) 


.5938(24) 
1.348(11) 


.5955(14) 
1.360(7) 


.6000(32) 
1.358(14) 


P 

L 


= 1.803: 

10 


20 


40 


80 


160 




9 
9r 


.5979(16) 
1.335(10) 


.5912(18) 
1.327(10) 


.5823(14) 
1.307(9) 


.5807(25) 
1.306(14) 


.5772(17) 
1.302(10) 





P = 1.805: 



L 


10 


20 


40 


80 


160 


9 
9r 


.5918(13) 
1.308(7) 


.5821(17) 
1.298(8) 


.5721(19) 
1.274(10) 


.5583(8) 
1.242(4) 


.5341(19) 
1.173(11) 



P = 1.82: 



L 


10 


20 


40 


80 


9 


.5566(8) 


.5196(8) 


.4578(28) 


.3533(42) 


9r 


1.213(4) 


1.119(4) 


.964(10) 


.652(8) 
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Tab. 2: 



Correlation length £, susceptibility x md renormalized coupling g^ in the 
high temperature phase of the icosahedron model 





1.470 


1.550 


1.610 


1.665 


1.707 


L 


80 


140 


250 


500 


910 


t 


11.203(5) 


19.627(12) 


33.655(21) 


63.628(33) 


122.09(16) 


X 


181.88(10) 


479.43(36) 


1228.04(91) 


3774.5(2.1) 


12009(17) 


9r(L) 


6.404(28) 


6.487(35) 


6.546(39) 


6.684(36) 


6.715(88) 


9r(°o) 


6.467(28) 


6.552(35) 


6.595(39) 


6.718(36) 


6.765(88) 


# runs 


200 


100 


100 


149 


28 



Tab.3: 

Correlation length £ and renormalized coupling gn in the high temperature 



phase of i 

(3 


Tie (03) mo 
1.5 


del (from ^ 
1.6 


1/ 

1.7 


1.8 


1.9 


1.95 


L 


80 


140 


250 


500 


910 


1230 




11.030(7) 


18.950(14) 


34.500(15) 


64.790(26) 


122.330(74) 


167.71(17) 


9r(L) 


6.553(16) 


6.612(15) 


6.665(14) 


6.691(15) 


6.737(21) 


6.792(40) 


9r(oo) 


6.616(16) 


6.668(15) 


6.730(14) 


6.733(15) 


6.792(21) 


6.853(40) 


# runs 


344 


370 


367 


382 


127 


68 
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